IMPROVED ASYMPTOTICS OF THE SPECTRAL GAP 
FOR THE MATHIEU OPERATOR 
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Abstract. The Mathieu operator 

L{y) ^ ~y" + 2a cos {2x)y, a G C, a =^ 0, 

considered with periodic or anti-periodic boundary conditions has, 
close to for large enough n, two periodic (if n is even) or anti- 
periodic (if n is odd) eigenvalues A~, A^. For fixed a, we show 
that 



A+ A- - ± 



■[(n-l)!p 



.4 



OO. 



This result extends the asymptotic formula of Harrell- Avron-Simon, 
by providing more asymptotic terms. 



1. Introduction 
The one-dimensional Schrodinger operator 

(1.1) L{y) = -y" + v{x)y, 

considered on R with vr-periodic real- valued potential v G ^^^^(R), is 
self-adjoint, and its spectrum has a gap-band structure (see Thm 2.3.1 
in [S], or Thm 2.1 in [IT]); namely, there are points 

A+ < < A+ < < A+ < A3 < A+ < A4 < A+ < ■ ■ ■ 

such that 

00 

Sp{L) = [j[K-i,K] 

n=l 

and the intervals of the spectrum are separated by the spectral gaps 

(-00, X^), (A^, A^), . . . , (A", A+), .... 

The points A~,A^ could be determined as eigenvalues of the Hill 
equation 

(1.2) - y" + v{x)y = \y, 

considered on [0,7r], respectively, with periodic boundary conditions 

(1.3) y{0) = y{7r), y'{0)=y'{7r), 

1 
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for even n, and with anti-periodic boundary conditions 
(1.4) y{0) = -y{7r), y'{0) = -y'{n), 

for odd n. See basics and details in P [171 flS] . 

The rate of decay of the size of the spectral gap 7„ = — is 
closely related to the smoothness of the potential v. Here we mention 
only Hochstadt's result [12] that an L^([0, 7r])-potential v is in C°° if 
and only if (7^) decays faster than any power of 1/n, and Trubowitz's 
result [25] that an L^([0, 7r])-potential v is analytic if and only if (7„) 
decays exponentially. 

If v is a complex-valued potential then the operator (11. ip is non- 
self-adjoint, so one cannot talk about spectral gaps. Moreover, the 
periodic and anti-periodic eigenvalues are well-defined for large n 
(see Lemma d] below) but the asymptotics of |A^ — A^| does not deter- 
mine the smoothness of f . In [22] V. Tkachenko suggested to consider 
also the Dirichlet b.v.p. y^n) = y{0) = 0. For large enough n, close to 
n"^ there is exactly one Dirichlet eigenvalue fin, so the deviation 

(1-5) 5n= |/i„-^(A+ + A;)| 

is well defined. Using an adequate parametrization of potentials in 
spectral terms similar to Marchenko-Ostrovskii's ones [181 US] ior self- 
adjoint operators, V. Tkachenko [221 121] (see also [23]) characterized 
C°°-smoothness and analyticity in terms of 6n and differences between 
critical values of Lyapunov functions and (—1)". See further references 
and later results in [21121 E] ■ 

In the case of specific potentials, like the Mathieu potential 

(1.6) v{x) = 2a cos 2a;, a 7^ 0, real, 
or more general trigonometric polynomials 

N 

(1.7) v{x) = 'y^^Ckexp{2ikx), Ck = c~k, < k < N < 00, 

-N 

one comes to two classes of questions: 

(i) Is the n-th spectral gap closed, i.e., 

(1.8) ^^ = A+-A- = 0, 

or, equivalently, is the multiplicity of A+ equal to 2? 

(ii) If 7„ 7^ 0, could we tell more about the size of this gap, or, for 
large enough n, what is the asymptotic behavior of 7„ = 'jn{v)7 
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E. L. Ince [H] proved that the Mathieu-Hill operator has only simple 
eigenvalues both for periodic and anti-periodic boundary conditions, 
i.e., 7n 7^ for every n e N. His proof is presented in [8]; see other 
proofs of this fact in [111 1201 IH] , and further references in [HI |26] . 

For fixed n and 0—7-0,0. Levy and J. Keller [16] gave an asymptotics 
of the spectral gap ■jn = 7n(a), v & (11. 6p : namely 

(1-9) 7n = A+-A; = |Mlll(i + 0(a)), a ^ 0. 



Almost 20 years later, E. Harrell [TO] found, up to a constant factor, 
the asymptotics of the spectral gaps of the Mathieu operator for fixed 
a as n — )■ oo. J. Avron and B. Simon [1] gave an alternative proof of E. 
Harrell's asymptotics and found the exact value of the constant factor, 
which led to the formula 

(1.10) -,„ = A:-A; = |M^(n-o(i)) n^oo. 

Later, another proof of (ll.lOp was given by H. Hochstadt [131. For 
general trigonometric polynomial potentials, A. Grigis [9] obtained a 
generic form of the main term in the gap asymptotics. 

In this paper we extend the result of Harrell- A vron-Simon and give 
the following more precise asymptotics of the size of spectral gap for the 
Mathieu operator (even in the case when the parameter a is a complex 
number) : 



;i.ii) A+-A, 



(a/4)" . a 



2 



n — > oo. 



Our approach is based on the methods developed in [H|5] , where the 
gap asymptotics of the Hill operator with two term potential of the 
form 

f(x) = Acos2x + 5cos4a;, A^O,Bj^O 

was found. The same methods and similar asymptotic estimates play 
a crucial role in the study of Riesz basis property of the root system of 
Hill operators with trigonometric polynomial potentials (see [7]). 

Let us mention that the paper [15] claims to provide "the formula 
which states the isolated terms of arbitrary number in the asymptotics 
of the sequence 7^." However, this claim is false due to a unavoidable 
technical mistake (in [15], formula (5) does not imply (4) for m = k + 1 
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since, by Stirling's formula, the remainder in (5) is much larger than 
the main term in (4)). 

To the best of our knowledge, f ll.lip is the first formula that gives 
more asymptotic terms than the formula of Harrell-Avron-Simon. 

2. Preliminaries 

Let Lper.+ (f) and Lp^r-iy) denote, respectively, the operator (11. ip 
considered with periodic (Per+) or antiperiodic {Per~) boundary con- 
ditions. Further we assume that v G L^([0,7r]) is a complex-valued 
potential such that 



(2.1) V{0) = / v{x)dx = 0. 

The following assertion is well-known (e.g., Proposition 1]). 

Lemma 1. The spectra of Lpf,^±{v) are discrete. There is an Nq = 
No{v) such that the union Un>NDn of the discs Dn = {z : |z — ?t,^| < 1} 
contains all but finitely many of the eigenvalues of Lpf,^.±. 

Moreover, for n > N the disc Dn contains exactly two ( counted with 
algebraic multiplicity) periodic (if n is even) or antiperiodic (if n is 
odd) eigenvalues A~,A+ (where Re X~ < Re or Re X~ = Re X^ and 
ImX- < ImX+). 

In view of Lemma [H 

(2.2) |A^-n2|<l, ioTn>No. 

Moreover, Lemma [T] allows us to apply the Lyapunov-Schmidt projec- 
tion method and reduce the eigenvalue equation Ly = Xy for A G Dn to 
an eigenvalue equation in the two-dimensional space = {L^y = n?y} 
(see [HI Section 2.2]). 
This leads to the following (see the formulas (2.24)-(2.30) in [6]). 

Lemma 2. In the above notations, A^ = + z, for \z\ < 1, is an 
eigenvalue of Lper±{v) if and only if z is a root of the equation 

(2-3) ^21 ^ _ ^22 = 0) 

where S^-^, S^'^, S"^^, S^'^ can be represented as 

oo 

(2.4) S'^{n,z) = Y,Sl\n,z), t,j = l,2, 
with 

(2.5) S',' = S',' = 0, S',' = Vi-2n), S',' = Vi2n), 
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and for each k = 1,2, 
(2.6) 

Sl\n,z)= Yl 
(2.7) 

Sr{n,z)= 

ji,...,jk¥'±n 

(2.8) 

Sr{n,z)= Y 

jl,...,jkJ^±n 

(2.9) 

Sl\n,z)= Y 

The above series converge absolutely and uniformly for \z\ < 1. 

Moreover, fl2.4p - fl2.9p imply the following (see Lemma 23 in [6J). 
Lemma 3. For any (complex-valued) potential v 
(2.10) S^\n,z) = S^\n,z). 





■n+ji)V{j2-ji) ■ 


■■VUk 


-jk-i)V{n 


-jk) 




( T)'^ — 9^ -1- 7\ ■ 




- j2 _|_ ^\ 




V{n 


+ Jl)^(j2-Jl)--' 




jk-i)V{-n 


- jk) 




[n? - jf + z) ■ 


■ ■ ■ (n2 - 


-fk + ^) 




v{- 


■n + ji)V{j2- jl)- 




-Jk-i)Vi- 


n - jk) 




(n2 - jl + z) 


■ ■ ■ {n^ 


-fk + z) 




V{n 


+ 3i)V{j2-3i)-- 




jk-i)V{n - 


Jk) 




(n2 - jI + z) ■ 


■ ■ (n^ — 


jI + 





Moreover, ifV{—m) = V{m) Vm, then 

(2.11) S^\n,z) = S^^{n,z), 
and if V{—m) = V{m) Vm, then 

(2.12) S^\n,z) = S^\n,z). 

Proof. For each /c G N, the change of summation indices ig = —jk+i-s, 
s = l,...,k proves that Sl^{n,z) = S'^{n,z). In view of (12. 4p and 
(1^ . d^TTU]) follows. 

In a similar way, we obtain that (12. lip and (I2.12p hold by using for 
each A; G N the change of indices is = jk+i-s, s = l,2,...,k. □ 

Further we consider only the Mathieu potential, i.e., 
(2.13) 

v{x) = 2a cos 2x = ae'^" + ae^", ^ (±2) = a, V{k) =0 if A; 7^ ±2. 

For convenience, we set 
(2.14) 

an{z) := S'\n,z) = S^'{n,z), /3n{z) := S^\n,z) = S'\n,z). 
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In these notations the basic equation (12. 3p becomes 

(2.15) {z - a,,{z)f = {(5^{z)f . 

By Lemmas [T] and [21 for large enough n G N, this equation has in the 
unit disc exactly the following two roots (counted with multiplicity): 

(2.16) z- = \--n\ z+ = \+-n\ 

3. Asymptotic Estimates for z^ and an{z). 

In this section we use the basic equation (I2.15P to derive asymptotic 
estimates for the deviations z'^. It turns out that |/3„(z)|, \z\ < 1, is 
much smaller than |a„(z)|, so it is enough to analyze the asymptotics 
of an{z^) in order to find asymptotic estimates for z^. 

The following inequality is well known (e.g., see Lemma 78 in [B]): 

(3,1) V_i--<ii£5^, forneN. 

Lemma 4. If \z\ < 1, then 



(3-2) V r-9 — ^ — -9 r < ( 

\^ 3i + z\- ■ ■ \n'^ — jf, + z\ \ 
Proof. If |z| < 1 and j ^ in, then 



4 log Qn 



n 



|2 '2i |\|2 '2| -i\^i2 "21 

\n — J + z\>\n — J | — l>-|n — j \ 



2' 



Therefore, 



V ^ < 2 

^-^ \'n? — ii + z\ - ■ ■ \'n? — il + z\ ~ 




ji,...,j„^±n 

SO follows from □ 

The next lemma gives a rough estimate for /3„(-2); we improve this 
estimate in the next section. 

Lemma 5. For \z\ < 1 we have 

logn 



(3.3) Pniz) = O 



n 



Proof. If z/ < n — 1, then all terms of the sum Sl^{n.,z) in (12. 9p van- 
ish. Indeed, each term of the sum Sl^{n,z) is a fraction which nom- 
inator has the form V {xi)V {X2) ■ ■ ■ V{xu^i) with xi = n + ji, X2 = 
j2 — ji, ■ ■ ■ , Xu+i = n — jy. Therefore, if z/ < n — 1 then there are no 
xi, a;2, . . . , Xy+i G {—2, 2} satisfying xi + X2 H h x^+i = 2n^ so every 
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term of the sum S'^^(n, z) vanishes due to fl2.13p . Hence, by f l2.13p we 
have 

oo 

\f^n{z)\< Yl E 

u=n—l ji,...,ji,^±n 

SO ([33D follows from □ 
Lemma 6. In the above notations, 





a 




\n'^ — ii + z 




?T,2 — j2 -|- 2; 



Proo/. In view of ([23D, ([MD and fl27[D . we have 

00 

(3.5) a„(z) = ^ylp(n,z), 

p=i 

where 
(3.6) 

V(-n + 3i)VU2 - 31) ■ ■ ■ VUp - 3p-i)V{n - jp) 



Ap{n,z)= 

ji,...,jp^±n 



(n^ - j( + 2;) ■ ■ ■ (n2 - j2 + z) 



First we show that 
(3.7) A2k{n,z) = V/ceN. 

Indeed, for p = 2k each term of the sum in (13. 6p is a fraction which 
nominator has the form V{xi)V{x2) ■ ■ ■ V{x2k+i) with 

Xi = -n + ji, X2 = 32 - 3i, X2k+i = n- 32k- 

Since Xi + X2 + ■ ■ ■ + X2k+i = 0, it follows that there is io such that 
XiQ 7^ ±2, so V"(xi„) = due to (I2.13p . Therefore, every term of the 
sum A2k{n,z) vanishes, hence (13. 7p holds. 

Next we estimate iteratively, in two steps, an{z) and z^. The first 
step provides rough estimates which we improve in the second step. 

Step 1. By (13. 6p . we have 

V{-n + 3i)V{n~3i) 



A(n,z)= Y 



- 3I + z 
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In view of f l2.13p , we get a non-zero term in the above sum if and only 
if ji = n + 2, or ji = n — 2. Therefore, 

(3.8) 

^2 ^2 

,2 



a? o 8 — 22; 



n2-(n-2)2 + 2^n2-(n + 2)2 + 2 (4^)2 _ (4 _ ^)2 ^ 
which imphes that 

(3.9) Ai{n,z) = O (^-^^ for \z\ < 1. 

On the other hand, from dSSD, (IS2D and dSSD it follows that 

(3.10) |A2._i(n,^)|<|ap^fi^V' \ A: = 2,3 



n 



which implies 

00 00 //I 1 « \ 2^-1 

(3.11) l^2^-(-' < E ) = o (l/n^) . 
Hence, by (13. 9p and (13. lip we obtain 

(3.12) an{z) = O [ij] for \z\ < 1. 



Furthermore, from fl2.15p . (12.160 and (13. 3p it follows immediately 
that 

' 1 



(3.13) ^±_«„(z±) = -), VfcGN. 



n 



Therefore, (I3.12p implies that 

(3.14) 4 = 
Step 2. By we have 

(3.15) A,{n,z) = ^^+o(-\ iiz = 0{l/n'). 



2n? 

Let us consider 

^ V- V{-n + 3i)VU2 - Ji)V{j3 - j2)V{n - J3) 
.i.W±n (n2-j2 + ;^)(n2-j| + z)(n2-j| + ^) • 

In view of (I2.13p . we get a non-zero term in the above sum if and only 
if 

ji = n + 2; j2 = n + 4; j3 = n + 2, 
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or 

ji = n - 2; 32 = n- 4; js = n - 2. 

Hence, 
A ( \ — 

^^"^'^^ ~ [n2 - (n + 2)2 + - (n + if + z][n^ - (n + 2)^ + z] 



^ [n2 - (n - 2)2 + z][n^ -{n- 4)2 + ^][n2 - (n - 2)2 + 

so it is easy to see that 

(3.16) A^in, z) = O (^-^^ if 1^1 < 1. 
On the other hand, by fl3.10l) we have 

oo oo //I 1 \ 2fc— 1 

(3.17) i^2.-i(-, -)i < E (^T^ ) = ^ (v-^) • 

fc=3 fc=3 ^ ^ 

Therefore, by fl37[5|) . fl3J6D and flSTTD imply that 

(3.18) a„(^) = ^ + (^^) if ^ = 0{l/n'). 
Hence, from fl3.13p it follows that 

(3.19) -^ = |^ + 0(i 



□ 



Remark. From (13. 8p and f l3.19p it follows that 

a2 a2 a'* _ / 1 



(3.20) ^,(„,.J) = --, + ^_-r_ + o^_^. 

Similarly, it is easily seen that 
(3.21, ^3(,,,J) = _|!_ + o(-i 

On the other hand, analyzing A^ln, z) one can show that 

(3.22) A^in, z) = O (J^^ if k| < 1. 
Moreover, by fl3.10p we have 

(3.23) f2\A2k-i{n,z)\=o(^\ if k| < 1. 

fc=4 ' 
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Hence, in view of f l3.13p . the estimates f l3.20p -( l3.23p lead to 

,2 ^2 / ]^ 



(3.24) 



a a 
2^^ 2^^ 



This analysis could be extended in order to obtain more asymptotic 
terms of z^, and even to explain that the corresponding asymptotic 
series along the powers of l/n contains only even nontrivial terms. 
However, in this paper we need only the estimate fl3.19p . 

The following assertion plays an essential role later. 
Lemma 7. With 7„ = — = — z~ 



(3.25) aniz-^) - an{z^ ) = in 

Proof. By (13. 5 p and (13. 7p we obtain 



,2 



in"- 



+ 0[- 



d 



dz 



a„{z) dz. 



(3.26) a„(4) - an{z-) = A^{n, z+) - A^{n, z^) + 
where 

din{z) = ttniz) - Ai{n, z) = Asin, z) + A^{n, z) A . 

In view of fl3A6D and fl3Tn) . 

Cin{z) = 0{l/n^) for \z\ < 1. 
By the Cauchy formula for derivatives, this estimate implies that 
dajdz = 0(l/n^) for \z\ < 1/2. 

Hence, we obtain 



d 



dz 



(3.27) 

On the other hand, by (13. 8p 

^i(^>2;+) - Ai(n,z, 



an{z) dz = 7„0 



1 



-2^^ 



-2z' 



(4^)2 _ (4 _ ^+)2 (4^)2 _ (4 _^-). 



In 



-32^2 - 32 + 8(2+ + z- 



[[{Anf-{A-ztY][{Anf~{A-z-f] 
Therefore, taking into account (13. 4p . we obtain 



(3.28) 



M{n,z1. ) - Ai(n,z„ 



In 



—a 



+ 



In view of ^M), the estimates (13:271) and fl3:28|) lead to ([325]). □ 
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4. Asymptotic formulas for Pri{zt) and 7^. 

In this section we find more precise asymptotics of /3„(z^). These 
asymptotics, combined with the results of the previous section, lead to 
an asymptotics for 7„. 

In view of (12.131) . each nonzero term in (12. 9p corresponds to a A;-tuple 
of indices (ji, jfc) with ji, . . . , Ja; 7^ ±^ such that 

(4.1) (n + ji) + (j2 - Ji) + ■ ■ • + (jfe - 3k-i) + {n- jk) = 2n 
and 

(4.2) n+jij2 - ji, ---Jk -jk~un-jk G {-2,2}. 

By (14.11) and (14.21) . there is one-to-one correspondence between the 
nonzero terms in (12.91) and the "admissible" walks x = {x{t))'l^^ on Z 
from —n to n with steps x{t) = ±2 and vertices jo = —n, jk+i = n, 

s 

(4.3) J, = -n + ^x(t) ^ in, s = l,...,k. 

t=i 

Let Xn{p), p = 0, 1, 2, . . . denote set of all such walks with p negative 
steps. It is easy to see that every walk x G Xn{p) has totally n + 2p 
steps because ^ x{t) = 2n. Therefore, every admissible walk has at 
least n steps. 

In view of ([23D, (S, fl^TT^ and ^HM . we have 
00 

(4.4) /3„(z) = ^(Tp(n, z) with ap{n, z) = ^ h{n,z), 
where, for x = {x{t))^^^, 

(4.5) h{x, z) — - - ^ ■ -- - t -, ^ TK— r 

with ji, . . . , jfc given by (113D- 

The set X„(0) has only one element, namely the walk 

(4.6) i = m)t=i, m = 2 vt. 

Therefore, 



(4.7) .„(„, .) = hK, r.) = 

with jk = —n + 2A;, = 1, ■ ■ ■ , n — 1. Moreover, since 

n— 1 

JJ (^2 _ ^ 2A;)2) = 4"-i[(n - 1)!]^, 
fc=i 
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the following holds. 

Lemma 8. In the above notations, 

4(a/4)" 



(4 



[(n-l)!]2- 



It is well known that 
(4.9) 



k=i ^ ^ 



n G N, 



where 7 is the Euler's constant. 
Lemma 9. In the above notations, 



a^logn a^7 



4^3 4^3 



+ 



(4.10) ao(n,4) = ao(n,0) 
Proof. By (14. 7p . we have 

(4.11) ao(n, 4) = ^o(n, 0) H ( ^ + ^2 _ + 



rt-l 



fc=i 



For simplicity, we set 



n'^-{-n + 2ky Ak{n - k) 



. Then, 



n-1 



'n-1 



log + bk)-' = - E + bk) = -J2^k + 0[J2\bk 



vfc=l 



fc=l 



Using (13.41) . we obtain 



n-l 



n— 1 



4A;(n - A;) 



k=l 



2^2 



ifc=i 



+ 



vfc=l 



2n I ^ A; 



2n2 



Ol4 



By dM]), it follows that 



k=l 



An^ 4^3 



On the other hand, by (137 



n-l 



fe=i 



^ [4A;(n - A;)]2 / 



IMPROVED ASYMPTOTICS 13 

Hence, 



which imphes f l4.10p . 



□ 



Next we study the ratio ai{n, z)/ao{n, z). 
Lemma 10. We have 

(4.12) ai{n,z)=ao{n,z)-^{n,z) 
where 

n-l 

(4.13) ^{n,z) =^ipkin,z) 

k=2 



with 

(4.14) ifkin.z 



a" 



[n2 - (-n + 2kf + z] [n2 - {-n + 2k - 2f + z\ ' 
Proof. From the definition of X„(l) and (14.41) it follows that 



n-l 



(4.15) ai{n,z)= ^ h{x, z) = ^h{xk, z), 

xeX„{l) k=2 

where Xk denotes the walk with {k + l)'th step equal to -2, i.e., 

Now, we figure out the connection between vertices of ^ and 
follows: 

jaixk) = < jk-liO a = k + l, 

[ja-2{0 k + 2<a<n = 2. 

Therefore, by 04. 5 p 

(4.16) hixk,z) = h{^,z) 



a" 



(n2-[j,_i(0F + z)(n2-[j,(0P + z)' 



Since jkH) = -n + 2k, k = 2, ... ,n - 1, (Hl5|l and fHlHj) imply 
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Lemma 11. In the above notations, if z = 0(1/?^^) then 

(4.17) $(n,z) = $(n,0) + (l/n^) 
and 

n-l 

(4.18) $*(n, z) := ^ |(^fc(n, z)\ = <l>(n, 0) + O (l/n^) . 

Moreover, 
(4.19) 

Proof. Since 



fc=2 



$(n,0) = ^ + ^^ + ^%^ + O(lK). 



4n3 



4n3 



^k{n, z) 



n -1 



1 + 



n -1 



1 + 



n2 - (-n + 2A;-2)2 



n2 - (-n + 2/c)2_ 
it is easily seen that 

ipk{n, z)/Mn, 0) = 1 + O (l/n^) if z = 0{l/n^). 
On the other hand, ipk{n,0) = 0(l/n2), so it follows that 
(^fc(n, z) - (^fc(n, 0) = 0) O (l/n^) = O {l/n') li z = O (l/n^) 

Therefore, we obtain that 

n-l 

5^ |(/.fc(n, z) - (/.fc(n, 0) I = O (IK) if z = 0(1/^2). 



k=2 



The latter sum dominates both |$(n, 2;) — $(n, 0)| and |$*(n, 2;) 
$(n,0)|. Hence, fliTTD and fOSD hold. 

Next we prove fl4.19p . Since 

n— 1 _2 

$(n,0) = 5^ 



k=2 



lQ{k - l)k{n - k){n + 1 - k) 



by using the identities 
1 1/1 1 



1 f 1 



1 



k{n — k) n \k n — k J ^ {k — l){n + 1 — k) n \k — 1 n + 1 — k 
we obtain 

.2 4 



(4.20) 



$(^,0) 



16^2 



i=l 
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where 



71—1 ^ n—1 



k=2 

n — 1 



k=2 



{n — k)(n + 1 — k) 



n — 1 -. n—1 



1 



(A; - l){n - k) 



k=2 ^ ' k=2 

The change of summation index m = n + 1 — k shows that D2{n) 
Di{n), and we have 

(4.21) D,(n) . g - i) . 1 - ^ . 1 - i + O (i) 

k=2 ^ ^ ^ ^ 



Moreover, since 



\fc=2 fc=2 / k=2 



by (14. 9 P we obtain that 



21ogn 2(7-1) 21ogn 



n 



n-2 



n 



(4.22) D3(n) 
Similarly, 

D,in) = J—ly^ + y — I — 

"" — 11 m ^-^ n — m — 1 



n-2 



n-2 



n 



\m=l 



n 



1 ^-^ m 



m=l 



and (14.91) leads to 
(4.23) 



21ogn 27 21ogri 
-D4(^) = \ \ ^ \- O 



n n 

Hence, in view of fOOjl -f iOS]) . we obtain flCTj) . 
Proposition 12. VFe /iave 



(4.24) 



/3„(z±) = ao(n,0) 



1 + 



O 



□ 



Proof. From ( ICTl . (HA2ll . ( liTTll and ( l419ll it follows immediately 
that 



cri(n, 2;^) + (Jo(n, z^) = ao{n, 0) 



^77,2 4n^ 



+ 



1 
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Since /3n{z) = Yl'^=o '^pi^^ view of f l4.10p to complete the proof it 
is enough to show that 

(4.25) f2 crp(n, z^) = ao{n, z^) O (-\ . 

p=2 ^ 

Next we prove fl4.25p . Recall that ap{n,z) = J2xeXn{p) ^i^^^)- Now 
we set 

(^*p{n,z)= \h{x,z)\. 
xex„{p) 

We are going to show that there is an absolute constant C > such 
that 

(4.26) a;{n,z^)<a;_,{n,z^)-^, peN, n > N^. 

Since ao{n,z) has one term only, we have aQ{n,z) = \ao{n,z)\. 

Let p G N. To every walk x G X„(p) we assign a pair {x,j), where 
X G Xn{p — 1) is the walk that we obtain after dropping the first cycle 
{+2, —2} from x, and j is the vertex of x where the first negative step 
of X is performed. In other words, we consider the map 

ip : Xn{p) — > Xn{p - 1) X J, / = {-n + 4, -n + 6, . . . , n - 2}, 

defined by v?(x) = (x,j), where 



x{t) 



x{t) if 1 < t < A; - 1 

x{t + 2) iik <t <n + 2p-2 



where k = min{t : x{t) = 2, x{t + 1) = —2} and j = —n + 2k. 
The map is clearly injective, and moreover, we have 

(4.27) hix, z) = hix, z)-— —— ; tt^ r. 

^ ' ^ ' ' ^ ' ' {n^ - p + z){n^ - {] -2f + z) 

Since the mapping Lp is injective, from fl4.14p . ( I4.18P and f l4.27p it follows 
that 

(4.28) a;{n,z)<a;_^{n,z)-<^\n,z). 

Hence, by <KW^ and (HT^ . we obtain that f lOB]) holds. 
From fl4.26p it follows (since al{n,z^) = \(rQ{n, z^)\) that 

(^*p{n,z^) < \aoin,z^) \ • ( ^ ) • 

Hence, (I4.25P holds, which completes the proof. □ 
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Theorem 13. The Mathieu operator 

L{y) = —y" + 2a cos {2x)y, a G C, a 7^ 0, 

considered with periodic or anti-periodic boundary conditions has, close 
to for large enough n, two periodic (if n is even) or anti-periodic (if 
n is odd) eigenvalues X~, X^. For fixed nonzero a G C, 



(a/4)" 



4^3 



+ 4 



n ^ 00. 



(4.29) A„ - ±1^^ _ 

Proof. The basic equation fl2.15p sphts into two equations 

(4.30) z - an{z) - Pn{z) = 0, 

(4.31) z - an{z) + /3n{z) = 0. 

In view of f l3.3p and f l3.12p . it foUows that for large enough n 

\z - an{z) ± (3n{z)\ < \z\ if \z\ = 1. 

Hence, for large enough n, each of the equations f l4.3Up and (I4.3ip has 
only one root in the unit disc due to Rouche's theorem. 

On the other hand, by Lemmas [1] and [21 for large enough n the basic 
equation has exactly two roots z~,z:^ in the unit disc, so either z~ is 
the root of fOOD and z+ is the root of (OT]) . or z+ is the root of fOOj) 
and z~ is the root of (I4.3ip . Therefore, we obtain 

Now, (13:251) and (^M> imply, with 7„ = A+ - A;;, 



7n 



1 + 



+ - 



±2ao{n,0) 



a 



Therefore, 

7„ = ±2(To(n,0) 

= ±2ao(n,0) 



1 + 



1 - 



a 



a 

An? 



O 



n- 

2 / 1 

+ o{- 



1 - 



O 



a 



An^ 



Hence, in view of (jM]), (l429ll holds. 



□ 
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